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We have developed a model for experiments in which the bias current applied to a Josephson 
junction is slowly increased from zero until the junction switches from its superconducting zero- 
voltage state, and the bias value at which this occurs is recorded. Repetition of such measurements 
yields experimentally determined probability distributions for the bias current at the moment of 
escape. Our model provides an explanation for available data on the temperature dependence of 
these escape peaks. When applied microwaves are included we observe an additional peak in the 
escape distributions and demonstrate that this peak matches experimental observations. The results 
suggest that experimentally observed switching distributions, with and without applied microwaves, 
can be understood within classical mechanics and may not exhibit phenomena that demand an 
exclusively quantum mechanical interpretation. 



I. INTRODUCTION 

In 1974, Fulton & Dunkleberger [1] demonstrated the 
way in which a biased Josephson junction could be ther- 
mally excited from its zero voltage state. More precisely, 
they conducted experiments in which the bias current 
was steadily increased from zero towards its critical cur- 
rent. In the absense of any noise, thermal or other- 
wise, the junction would not switch until the bias cur- 
rent reached the critical value. However, with thermal 
noise, junctions were observed to switch with high prob- 
ability at bias currents that were slightly less than the 
critical value. This type of experiment has proven to 
be an extremely useful tool for probing the details of ef- 
fective models of the junctions themselves. Later work 
by Voss & Webb [2J extended the experiments to much 
lower temperatures and found what was interpreted to be 
evidence that the junction had entered a mode where an 
escape might be treated as Macroscopic Quantum Tun- 
neling (MQT) out of the effective potential well associ- 
ated with junction phase dynamics. 

Once the idea of macroscopic quantum behavior for 
Josephson junctions at low temperatures became ac- 
cepted, the possibility of the manifestation of discrete 
quantum levels within the effective potential wells re- 
ceived attention. The first experiment to consider this 
propostion was reported by Martinis, Devoret and Clarke 
(MDC) [3j. In that experiment, microwaves were di- 
rected onto a junction and the bias current was swept as 
before. The idea was that the microwaves would excite 
transitions from a lower to a higher level within the well, 
and that the macroscopic quantum variable - the junc- 
tion phase - would then tunnel out of the well from that 
higher level, resulting in an escape from the zero voltage 



state. Single and multiple peaks in the escape distribu- 
tions of two samples were reported in their experiment. 
These data were interpreted as signatures of the antici- 
pated level transitions dictated by quantum theory. 

Previoisly we analyzed the issue of the classical reso- 
nant frequencies in wells under both harmonic and an- 
harmonic approximations and found [4| that the classical 
theory gave results in good agreement with the experi- 
ments of MDC, suggesting that the classical Resistive and 
Capacitive Shunted Junction (RCSJ) model for a Joseph- 
son junction might have been dismissed prematurely in 
favor of the macroscopic quantum picture [5|. The im- 
portance of the anharmonic component to the potential 
well is evident from the fact that a measurement of es- 
cape cannot be realized in a harmonic well and, thus, one 
should not expect meaningful agreement between switch- 
ing experiments and classical analysis of plasma oscilla- 
tions. This notion was first presented in Ref. [6|, where 
the anharmonic theory successfully compared to accom- 
panying experiments on both direct and harmonic reso- 
nant switching. It should be mentioned that the anhar- 
monic classical RCSJ approach has since also produced 
good agreement with other experimentally observed fea- 
tures, such as Rabi oscillations and Ramsey fringes [3, [8j 
as well as tomographic reconstruction of anticipated den- 
sity matrices for a pair of capacitively coupled Josephson 
junctions [9j], which originally had been interpreted ex- 
clusively in terms of quantum entanglement. 

While the non-quantized RCSJ model has proven to 
replicate the primary resonant features of the experi- 
ments, available experimental reports on switching dur- 
ing bias-sweep contain details not yet directly analyzed 
using the non-quantized approach. One is the saturation 
of the width of the switching distribution as the ther- 



modynamic temperature is lowered, interpreted as a sig- 
nature of a quantum crossover temperature [2]; another 
is the set of multiple of resonant switching peaks, inter- 
preted as signatures of quantized energy transitions in 
the potential well |3j- 

In this paper, we reconsider the evidence for the macro- 
scopic quantum tunneling interpretation of these experi- 
ments. Inparticular, in the spirit of Kramers' statistical 
analysis (lOj, we develop a simple model of the swept- 
bias type of experiment based on classical thermal acti- 
vation, and show that it gives an excellent accounting of 
the peaks observed in a number of key experiments. 



II. MODELING A JOSEPHSON SWEPT-BIAS 
EXPERIMENT 

We consider Josephson junctions to be characterized 
by a supercurrent Iq sin cp with a critical current Iq and 
junction phase ip, and a capacitance C. The ass ociated 
junction plasma frequency is uoj = ^/2irlc /&oC, where 
$o = h/2e is the flux quantum. We assume junctions 
whose physical dimensions are much smaller than the 
Josephson penetration depth. 

From this perspective, the phase dynamics of a single 
Josephson junction subjected to a dc bias is equivalent 
to the motion of a particle on a washboard potential [U( . 
The particle sits in a potential well which becomes shal- 
lower at high bias currents. This in turn leads to an 
enhanced probability of noise activated escape. If the 
bias current is swept from towards Ic, then at some 
moment the junction will be observed to escape from its 
potential well and switch to a running state with finite 
voltage. In the experiments to which we will refer, this 
process was repeated many times in order to acquire a 
statistical profile of the distribution of bias currents for 
which the escapes from the zero-voltage state occur. 

For our numerical simulations, we imagine an equiv- 
alent scenario. Suppose there is an ensemble of M 
Josephson junctions. The bias on all junctions starts at 
and is incremented in N steps, with each step of du- 
ration At = (Nfs) -1 where fs is the sweep frequency. 
Each step is assigned a channel and the total counts in 
that channel indicate how many junctions have switched 
to a finite voltage state (escape from the potential well) 
during that interval. As the bias sweep proceeds, the 
original ensemble will have lost e\ junctions in the first 
interval, e2 junctions in the second interval, and so forth. 
Consequently, at the beginning of the n th bias interval, 

n-l 

there will be M — 2_, e j junctions not yet escaped. The 

i=i 
number from this remaining pool of junctions that will 

escape during the next At seconds will be 



M~J2 e i 



T(t n )At, n = 2,3,..JV (1) 



where T(t n ) is the probability of escape per unit time 
in the n th interval, Of course, the initial interval just 

satisfies 



ei = MT(t!)At 



(2) 



Equations (jl!2p will mimic a swept-bias experiment pro- 
vided a suitable expression is available for the escape rate 

r. 

In Kramers' theory [10], the thermal escape rate can 
be expressed 



r(t n ) = / n exp 



AU n 
k B T 



(3) 



where f n is the plasma frequency for the well specific to 
the n th bias interval and AUn/ksT is the height of the 
potential barrier divided by the mean thermal energy. 
Voss and Webb [2[ assumed an escape rate in this form. 
However, there has been an ongoing discussion regard- 
ing the suitability of the Kramers expression, and as De- 
voret, Esteve, Martinis, Cleland, and Clarke [12], and 
others, have pointed out, a better equation for the es- 
cape rate, from Biittiker, Harris, and Landauer (BHL) 
[I|, is 



r B HL(t n ) = a t / n exp 



' kT 



where 



a t 



\a 



i + V( 1 + ^€) 



(4) 



(5) 



According to Devoret et a/.|l2[, a = 1.4738. In this ex- 
pression, Q is a parameter that quantifies the dissipation 
in the junction; the lower the dissipation, the larger the 
Q. Devoret et at. [14] noted: "The prefactor depends 
only weakly on Q" and, they consequently used an es- 
cape rate in the form of Eq.(3). Similarly, Devoret et al 
in [12| stated that the value of the prefactor a t is "close 
to unity". 

For these reasons, we proceed with our classical sim- 
ulations of swept bias experiments using expression (|3j) 
for the escape rate. We shall consider the second-order 
effects of Q in Section VI. 

In the harmonic approximation 



fn = fj Vl^n 



(6) 



with r/ n being the normalized bias current within the n th 
bias interval. 

Combining these expressions, we obtain 



T(t n )At 



fj 

Nfs 



V^7n 



exp 



AU n \ 
k B Tj 



(7) 



which is required in Eq.(pQ). The prefactor before the 
exponential represents the number of plasma oscillations 



that can fit within the time window of the n th data ac- 
quisition channel - that is, the number of attempts that 
will occur in that time window. 

According to MQT theory at sufficiently low temper- 
atures tunneling of the phase variable becomes the dom- 
inant escape mechanism, in which case the escape rate 
takes the form 



where 



a q f n exp 



a a = 



-7.2 



MJn 
hfn 



0.87 



1207T 



7.2AU n 

hfn 



(8) 



(9) 



Clearly, in contrast to the classical expression Eq.(j3]), this 
escape rate does not depend on temperature. Voss and 
Webb [2] pointed out that the transition from classical to 
quantum behavior should take place around a " crossover 
temperature" satisfying 



7k B 



(10) 
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FIG. 1: (color online) Simulation results for a swept-bias 
experiment. The ten classical peaks had f3 values (left to 
right): 0.0551, 0.0375, 0.0259, 0.0169, 0.0123, 0.0092, 0.00656, 
0.00372, 0.00223, 0.00160. The single escape peak labelled 
MQT was computed using expression ((8} for the escape rate 
with Q = 50. 



III. MICROWAVES OFF 

The height of the barrier at the n th step is given by 
the well known expression 



AU n = 2 



2tt 



V 7 ! - r]l - T] n COS 1 T] n 



(11) 



so 



AU n 2 r 



k B T 



y/l-rfo- Vn cos 1 rj n (12) 



where 



>-(^)£ 



1. Voss & Webb (1981) 

As an example of an experimental simulation, param- 
eters were chosen to be: N = 5000, M = 100,000, 
fj = 35.53 GHz, fs = 10 Hz, numbers consistent with 
the experiments of Voss and Webb [2| (although they 
gave no value for the number of channels in their system). 
The evolution of the peaks shown in FigJT]is in very good 
agreement with the experimental data in Fig.l of Voss 
and Webb [2]. Also shown in FigJT] is the single MQT 
peak from a simulation using the escape rate expression 
Eq.flSJ) with Q = 50. The conclusion to be drawn is that 
when the value of j3 drops below a crossover equivalent, 
macroscopic quantum behavior should take over from the 
classical escape process and the temperature independent 



peak marked MQT should become frozen in place. In 
such a case not only the peak widths but also the peak 
positions must remain constant. Then, none of the clas- 
sical peaks to the right of MQT would be observed in an 
experiment. 

The crossover temperature, Eq.JTDJ), is a function of 
the natural frequency of a particular well, and this in 
turn is controlled by the applied bias current as speci- 
fied in Eq.(j6]). Therefore, in a swept-bias experiment 
one is also sweeping the natural frequency of the contin- 
uously varying well shape. It is simple to use Eqs. (|6I10|) 
to plot the dependence of crossover temperature on bias 
current; this is shown in FigI2J The positions (bias val- 
ues) of the experimental peaks, indicated by solid dots, 
were manually extracted from Fig.l in [2[ using digitizing 
software [15]. A vertical line marks the point at which 
the sample temperature has dropped below the crossover 
characteristic - this occurs at T w 155 mK and a bias 
of 1.470 jiA. But note that the two experimental escape 
peaks for temperatures T = 95 mK and T = 5 mK are 
below the anticipated quantum transition temperature 
and appear not to have frozen at 1.470 \iA, but instead 
continue to advance beyond the MQT stopping point, 
into the shaded 'forbidden' zone. While this progression 
of escape peaks towards higher bias values is contrary to 
the expectations of the MQT model, it is consistent with 
the classical escape model. 

The apparent saturation of the widths of the escape 
peaks below the crossover tempereature, noted in [2] was 
claimed to constitute the "first compelling evidence for 
the existence of quantum tunneling of a macroscopic vari- 
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the observed temperature dependence of the peak widths 
in [2[ could be replicated using T/v = 63 mK and with 
classical escape theory alone. Such an elevated sample 
temperature, possibly due to self-heating, was also noted 
in 03 

We have run our simulation with T replaced by T e ff 
and T/v = 63 mK and the results, virtually identical to 
those in [16[, are shown in Figj3] By most standards, the 
agreement between experiment and classical theory is ex- 
cellent. A vote in favor of a quantum signature in these 
data could only be supported by a much finer analysis 
including error bars in the experimental peak widths. As 
it stands, the classical model fits the experimental data 
exceptionally well over two orders of magnitude in bath 
temperature with the above-mentioned suggested effec- 
tive sample temperature. 



FIG. 2: (color online) Dependence of the crossover tempera- 
ture on bias current, for sample parameters of Voss and Webb 
Q] • Dots mark the positions of the seven lowest temperature 
peaks shown in Fig.l of [2]. At 1.470 fiA the experimental 
data drop below the crossover boundary. According to the 
MQT hypothesis, the escape peaks should then become tem- 
perature independent and peak positions would be expected 
to lie on the vertical line, as depicted by open circles. 
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FIG. 3: (color online) Dots: experimental peak widths from 
Fig. 3 of ref. 2]. Dashed line: classical simulation with an 
effective sample temperature of T + 63 mK where T is the 
bath temperature. 



able". Using digitizing software, the experimental data 
oints for the peak widths were extracted from Fig. 3 in 
. These points are plotted in Figj3] More than twenty 
years ago, Cristiano and Silvestrini [16| proposed that the 
presence of some additional noise could raise the sample 
temperature above the bath temperature T such that 
T e ff = T + Tat. In particular, they demonstrated that 



IV. MICROWAVES ON 

We now consider appropriate modifications to the 
swept-bias model needed to take into account the action 
of applied microwaves. In such experiments, a fixed fre- 
quency microwave source is used and the bias current is 
ramped up. 

Our numerical solutions of the equation of motion of 
the phase [ll| with both dc and ac bias reveal induced 
oscillations around the minimum in the potential well. 
The amplitude of these oscillations depends on the dc 
bias chosen. There will be a particular bias, r] res , at 
which the the oscillation amplitude is a maximum. The 
question is: What is the value for r] res ? 

In the harmonic approximation, the natural frequency 
of a well is related to the dc bias through Eq.((6j). How- 
ever as shown in [4, 6], when the amplitude of the phase 
oscillations is large, the cubic nature of the well comes 
into play and an anharmonic approximation takes the 
place of Eq.(j6]); it is. 



In — J 3 



\ 



(Jo(A) + J 2 (A)) 



Vn V 



MA) J 



(14) 



where J p is the Bessel function of pth order, first kind 
and A is the amplitude of the oscillation. Thus, the 
resonance frequency is depressed for increasing oscilla- 
tion amplitudes A. It has been found that situations, 
in which resonant states produce non-zero and non-unity 
switching probabilities, are given for oscillation ampli- 
tudes near the inflection point of the potential well. In 
the limit rj — > 1, this value of A is given by the explicit 
expression A 2 « ~(1 — r] res ). 

Setting f n = f ac in Eq.(|6j) would give one answer 
for the resonant bias r] resi while Eq. (fT3|) would give a 
slightly smaller answer. This means that without knowl- 
edge of the strength of the microwaves at the junction, 
the best one can say is that r] res must lie somewhere be- 
low the value for A = (as also seen experimentally in 



|18| ) and in the vicinity of the interval spanned by the 
two values for A = and A 2 « |(1 — r] res ). As an exam- 
ple, for a microwave frequency f ac /fj = 0.350, the limits 
of r] res are: 0.9925 and 0.9887. 

Let the amplitude of the induced phase oscillation to 
the right of the minimum point of any well be denoted 
Scp n . For not too large excitations, 5(p has a bell-shaped 
distribution, centered at r] res . At r] res the ac field is trans- 
ferring a maximum amount of energy into the junction. 
On either side of this optimum bias, the amplitude of the 
phase oscillations diminishes and the absorbed energy de- 
clines. For a bell shaped distribution of Sep we used the 
following heuristic expression, 



fyn 



(Vn ~ Vres) + b 2 



(15) 



There are two parameters here: b which sets the sharp- 
ness of the distribution, and a which sets the peak value. 
Whenever sustained phase oscillations are induced, the 
added energy is 



Mhn 

k R T 



= P 1 {[- cos </? max - ?7„</? max )] 

- [- COS(p m i n - Tin (<£min)]} 
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J l.l l i 



(16) 




with cp min = sin : r\ n and (p max = (p min + 5<p n . 
This will reduce the escape barrier in the n th interval to 
an effective value 



AUeff _ A/7 n AU ln 
k B T ~ k B T k B T 



(17) 



and this is what thermal noise needs to overcome. This 
effective barrier height replaces the original in Eq. ([7]) and 
then the simulation can proceed as before. 

We now apply this simulation to several sets of pub- 
lished data. First we consider experiments which showed 
only a single microwave induced escape peak. 



A. Single Microwave Induced Peaks 

As noted already, swept bias experiments yield his- 
tograms for the escape probability. Some authors prefer 
to convert such data to escape rates T as a function of 
bias current. Then a relative rate, with and without mi- 
crowaves, [r(P) - r(0)] /r(0) may be plotted. This has 
the effect of stripping away the thermal escape peak (as 
discussed in Sections 3 & 4), thereby isolating purely mi- 
crowave induced phenomena. 



1. Martinis, Devoret & Clarke (1985) 

Consider the results presented in Martinis, Devoret, 
and Clarke [3[. Their Fig. 3 is reproduced in the up- 
per panel of FigHJ The junction was characterized 



FIG. 4: (color online) Comparison of the experimental results 
of Martinis, Devoret, and Clarke [3] (upper) with our simula- 
tion (lower). 



by the following parameter values: Iq = 9.489 jiA 
and C = 6.35 pF. This gives a junction zero bias 
plasma frequency fj = 10.72 GHz; hence the microwave 
frequencies of 3.7, 3.6, 3. 5, and 3AGHz correspond to 
0.3451,0.3358,0.3265 and 0.3172 in dimensionless form. 
This experiment was carried out at T = 18mK. 

For comparison of these experiments with classical re- 
sults, we simply make use of the escape rate expression 
Eq.([3]) with a barrier given by Eq. ([T2]) in the absense 
of microwaves, or with a reduced effective barrier given 
by Eq.([T7]) when microwaves are present. For each of 
the four microwave frequencies, the anharmonic result 
Eq.(p3j) was used to obtain the bias current r] res which 
selects the well that is resonant with the excitation. With 
this normalized bias, Eq. (fT6|) together with Eq. (fT5|) per- 
mits the reduced barrier height to be calculated. 

The remainder of the parameters chosen to match the 
situation were: a = 0.10, b = 0.0010 and /? = 0.0000796. 
The results for the classical escape rate calculations are 
shown in the lower panel of Fig 21 As can be seen, the clas- 
sical results agree very well with the experimentally ob- 
served peaks in the escape rates for these four microwave 
frequencies. 
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FIG. 5: (color online) Comparison of experimental data from 
Thrailkill et al. 19] (upper) and our simulation results (lower) 
where closed circles indicate results without microwaves. 



2. Thrailkill, Lambert, Carabello & Ramos (2009) 



B. Multiple Microwave Induced Peaks 



Every swept bias experiment with microwaves present 
yields at least one ac induced escape peak. This includes 
both Figs.2 and 3 in [3], Fig.3 in [17], Fig.2 in [19], and 
Fig. 6. 5 in [20]. But only in two instances was a second 
microwave peak in evidence. 

In Fig.2, of [3], one of the experimental samples exhib- 
ited an additional microwave escape peak. For that ex- 
periment, the sample parameters were: Iq = 30.572/iA, 
C = 47pF and so fj = 7.07GHz. Hence the normalized 
microwave frequency of2GHz was 0.2829. The two peaks 
were at 30.43424/M and 30.41391/^4, which in normal- 
ized units are 0.995490 and 0.994829. These two points 
are included in Figj6j they are both quite close to the 
anharmonic curve. 

In Fig.3 of |17j, there are two microwave induced es- 
cape peaks. The sample parameters were: Iq = 14.12 \±A, 
C = 4.2 pF yielding a junction plasma frequency fj = 
16.1 GHz. The microwave frequency was 5.7 GHz which 
is 0.354 in normalized units. The two peaks are at bias 
currents of / = 13.9907/iA and / = 13.9530/iA, which 
in normalized units are 0.9908 and 0.9882. These two 
points are included in FigHJ Interestingly, one appears 
to be on the anharmonic curve, while the other is close to 
the harmonic curve, perhaps suggesting that the system 
can resonantly respond to either condition. The pos- 
sibility for multiple states in the ac driven anharmonic 
potential is consistent with previously reported observa- 
tions (see Fig.l in Ref.[8]). It should be pointed out 
that the solid line Berkley's Fig.3 was described as "a 
Lorentzian fit to two peaks", meaning it is not in any 
sense a theoretical prediction and so does not constitute 
confirmation of quantum expectations. The labelling of 
the two peaks as |1) — » |2) and |0) — » |1) is based on 
assumptions regarding applicable physics. 



Thrailkill et al. [19| carried out swept-bias experiments 
on a Josephson junction characterized by Iq = 9.485 [iA 
and C = 4.7 pF. The Josephson plasma frequency was 
thus fj = 12.46 GHz. The escape probability distribu- 
tions were measured at several different combinations of 
temperature and microwave frequency, as shown in the 
upper panel of Figj5] For the classical simulations, the 
parameter values were: a = 0.09, b = 0.003 (equivalent 
to a half width in the distribution of 0.028 fiA) and from 
top to bottom, r] res = 0.972, 0.976, 0.980, 0.984. The sim- 
ulation results clearly are in very good agreement with 
the experimental data. Note that the ac resonance and 
the thermal peak are nearly on top of each other, so the 
peaks almost merge. Also, the simulation exhibits the 
same shifting effect as in the experiments - the thermal 
peak without microwaves becomes displaced slightly to- 
wards lower bias values when the microwaves are turned 
on. 



In the previous section, it was noted that the four sin- 
gle peaks, each at a different microwave frequency, from 
Fig.3 of Q are reasonably reproduced by classical calcu- 
lations of relative escape rates. These classical points are 
shown in Figj6]as diamonds that lie along the anharmonic 
curve. However, the comparison of experiment with clas- 
sical and quantum theories described in [3[ contained the 
statement: "Furthermore, the measured positions of the 
resonances are clearly very different from a classical pre- 
diction for the resonant activation of the particleoscillat- 
ing at the plasma frequency (dashed line)" . This dashed 
line appears in part (b) of the figure and is in fact the 
harmonic approximation. But we see from our Figj6]that 
the fair test of the classical model is the anharmonic ap- 
proximation, and at the very least the classical model is 
as successful as the quantum hypothesis. 
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FIG. 6: (color online) Comparison of classical model and 
quantum model with experimental data. Circles are exper- 
imental points, squares are predictions of the quantum model 
(digitized from the positions of the arrow markers in Fig. 3(b) 
in [3]), and diamonds are predictions of the classical model us- 
ing the anharmonic approximation for resonance frequencies 
Eq. ([14[) . The harmonic approximation, Eq.JB} is also plotted. 
The bottom pair of points are taken from the experimental 
data in Fig. 2 of [3]; the upper pair of points (stars) are the 
two observed peaks from Fig. 3 in [171 ]. 



V. EFFECT OF SWEEP FREQUENCY 

It is important to note that peaks in the escape prob- 
ability distributions are not like lines in atomic spec- 
tra in that they are a manifestation of both the fun- 
damental physics associated with escape rates and the 
way the experiment is performed - specifically, the fre- 
quency at which the bias current is swept from zero to 
its critical value, (fs)- This issue was addressed in a 
swept bias simulation with: N = 50, 000, M = 100, 000, 
fj = 7.072 GHz, f s = 10 Hz and /? = 0.00129. The 
results for microwaves OFF are shown in FigH Clearly, 
the exact location and shape of the thermal escape peak 
are determined in part by the speed with which the bias 
current is ramped towards the critical value. As might 
be anticipated, increasing the sweep rate moves the peak 
towards higher bias currents. 



VI. EFFECT OF DISSIPATION 

To illuminate some aspects of the effects of dissipation, 
we carried out swept bias simulations using Eqs. (j4]) and 
([5]), and for this example with parameters Iq = 14.12/iA, 
C = A.2pF, f3 = 0.00018. The microwave frequency for 
the was set at b.7GHz and the assumed excitation pa- 
rameters were a = 0.085, b = 0.0012. Typical values of 
the dissipation Q for underdamped Josephson junctions 
lie in the range 20 to 50. Escape histograms were re- 
peated for a number of choices of the dissipation constant 
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FIG. 7: (color online) Simulation results showing the effects 
of bias sweep frequency fs on the position of the peak in the 
escape probability distribution. 
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FIG. 8: (color online) Swept bias simulations with selected 
values of the parameter Q u sing the escape rate due to 
Biittiker, Harris, and Landauer[13(]. 



parameter and the results are shown in Fig.??. Note that 
the microwave peak does not shift, but the position of the 
purely thermal peak varies with Q. Therefore, the partic- 
ular value of the junction dissipation might have a slight 
effect on predictions of the Voss & Webb type of experi- 
ment, but would not influence the positions of microwave 
induced peaks. 



VII. DISCUSSION 

The decades old papers of Voss and Webb, and Marti- 
nis et al. appeared to convincingly demonstate the (antic- 



ipated) appearance of MQT in superconducting circuits 
operating below a crossover temperature. The classical 
model was subsequently discarded as a possible source for 
observed phenomena at millikelvin temperatures. In this 
paper we have shown that this assertion may not be jus- 
tified, and that these early foundational experiments can 
certainly be modeled successfully within a purely classi- 
cal device description. 

With respect to the experiments of Voss and Webb, we 
have demonstrated that there was not strong evidence 
that the junction had entered a macroscopic quantum 
state even at the lowest temperatures. A classical model 
with some self-heating gives a more consistent description 
of those observations. 

With respect to the experiments of Martinis et a/., we 
have demonstrated that in the presence of microwave ir- 
radiation the additional peaks which appear in the swept 
bias escape distributions are just as well accounted for 
within the classical resonant activation model as by the 



proposed macroscopic quantum model. 

The key issue in this situation was nicely expressed by 
Devoret, Martinis, and Clarke [2l| as follows: "An ex- 
periment cannot prove a theory, but only invalidate an 
alternative theory." . The present study should therefore 
be seen in this context - the classical theory for these sys- 
tems has not yet been ruled out. Therefore, an exclusive 
presumption of MQT in these systems is not justified. 
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